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1. Introduction 

Quantum graphs have become a widely employed model for a range of phenomena in 
mathematical physics from Anderson localization and quantum chaos to problems in 
mesoscopic physics, microelectronics and the theory of waveguides, see e.g. [271 EHj- 
Much of the success for the analysis of these problems on graphs can be attributed to 
the existence of well behaved spectral functions. Consequently spectral determinants 
and trace formulae in graphs have become important problems widely studied in the 
literature [2l|3liai3[71[9l[8l[T0l[IIl[l3l[l5l[T8l[T9ll^ Here we consider 

the Schodinger operator —A + V{x) on a metric graph where each bond of the B 
bonds is associated with an interval [0,Lfe], b = 1, ... ,5. Matching conditions at the 
vertices of the graph are parameterized by 2B x 2B matrices A and B according to 
the scheme introduced by Kostrykin and Schrader [21] , see Section [21 We also restrict 
the possible matching conditions to local matching conditions which both respect the 
graph topology, only relating values of the functions and their derivatives at the ends 
of intervals where they meet at a vertex, and are independent of the metric structure of 
the graph, namely the bond lengths {Lb}b=i,...,B and the energy. The matrices A and B 
can thus be considered to define the graph topology as they record how the bonds are 
connected. 

In particular we will be concerned with the recent broad conjecture, introduced by 
one of the authors in [32], for the spectral determinant which is defined formally as 

oo 

^(7) = det[7 - A + V{x)] = JJIt + Ej) , (1) 

3=0 

where < Eq ^ Ei ^ . . . is the spectrum of the Schrodinger operator and 7 
is some spectral parameter; the condition on the positivity of the spectrum can be 
relaxed by accepting additional technicalities [21]. The first result for the spectral 
determinant of the Schrodinger operator on a graph with general vertex matching 
conditions was obtained by Desbois [9]; the method was, however, unable to determine 
the 7-independent prefactor of the determinant. 

In order to regularize the formal definition of iS'(7) we employ the spectral zeta 
function, 

00 

as) = Y,ii + ■ (2) 

i=o 

The zeta function of the Schrodinger operator can be formulated by extending a 
technique introduced for graph Laplacians by some of the authors [T71 [18]. Having 
derived the zeta function the regularized spectral determinant is then defined by 
= 6xp(— C'(0)). More precisely we prove the following theorem along the lines 
of the conjecture. 

Theorem 1. For the Schrodinger operator on a graph with local vertex matching 
conditions defined by a pair of matrices A and B, with AB"!" = BA"!" and rank(A, B) = 2B, 
and potential functions Vb{xb) G C°° for b = 1, . . . ,B the spectral determinant is 
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\b=i 




-2 



) 



det (A + BM(7)) 
Cat 7^ 



fbi^'i 7) ^(7) cire defined in terms of tlie solution to a boundary value problem 
on the interval [0, Lh] described in Section HI is the coefficient of the leading 
order t — )■ 00 asymptotic behavior of det (A + BM(t^)), determined in Section l5.ll 
Furthermore the exponent P characterizes the t — ?■ behavior of the determinant; 
P = if det (A + BM(0)) ^ 0, the generic case, otherwise det (A + BM(7)) ~ 'y^ as 
7^0. 

To complete the structure of the article, in Section [2] we define the terminology of 
the quantum graph model. Section [3] introduces the technique, followed subsequently, 
in the simple case of the Schrodinger operator on an interval with Dirichlet boundary 
conditions and a graph with Dirichlet conditions at the ends of every bond. To derive 
the zeta function of a general graph we start from a secular equation whose solutions 
are ^/E^', an appropriate secular equation for the Schrodinger operator is identified in 
Section m In Section we formulate the graph zeta function with general local matching 
conditions at the vertices. Finally, in Section |6l we apply the results for the zeta function 
to prove Theorem [TJ The Appendix shows the derivation of the t — ?■ 00 asymptotics of 
the function 

2. Graph model 

A graph G is a collection of vertices v = 1, . . . ,V and bonds b = 1, . . . ,B. Each bond 
connects a pair of vertices b = (f , w) and we denote by o{b) = v the initial vertex of b and 
t{b) = w the terminal vertex, b = {w,v) will denote the reversed bond with the initial 
and terminal vertices exchanged. We consider undirected graphs with B undirected 
bonds, so b and b refer to the same physical bond but the availability of two labels for 
each bond will be notationally convenient. To be explicit the label b will refer to the 
bond with o{b) < t{b). m„ will denote the number of bonds meeting at a vertex v, the 
valency oi v. To determine the valency b and b refer to a single bond. 

In order to define a metric graph we associate with each bond b an interval [0, L^] 
where is the length of b and xj, = at 0(6) and Xh = Lb at t{b). So Xb is the distance to 
a point in the interval measured from the initial vertex o{b) and we also use = Lb — Xb 
for the distance to the same point measured from the terminal vertex of b. The total 
length of G is denoted by C = X]fc=i ^b- A function on G is defined by specifying 
the set of functions {ipb{xb)}b=i,...,B on the collection of intervals. The redundancy of 
notation enforces the relation ipbi^b) = ipbi^b) — i^i^b — ^b)- The Hilbert space of G is 
consequently 



Motivated by physical applications we consider here Schrodinger operators on the 



B 



'H = ^L\%Lb]). 



(3) 



6=1 



Spectral determinants on graphs 



4 



metric graph. The eigenproblem on bond b is 



d 



^doT 7 Vb{xb)ilJb{xb) = k ipbixb) ■ (4) 

The set {Ai, . . . ^Ab} defines a vector potential on the graph constant on each bondj^ 
For consistency Aj^ = —Ab as the direction is reversed when changing coordinate. 



3. Dirichlet determinants 



3.1. Laplace operator on a finite interval 

We first analyse the determinant of the Laplace operator on a single finite interval [0, L] 
with Dirichlet boundary conditions, a wire. The spectrum of the wire is Ej = fo'^ 
j G N* and we can compute the zeta-function directly, 

oo / J \2s 

CDi.(.) = 5^E7«= - W^S). (5) 

j=l ^ ^ 

Differentiating, 

C^-M = 2Cj,{0) + 2Ci?(0) log(L/7r) = - log 2 - log L , (6) 

where we have used C-r(O) = — | and Ci?(0) = — |log(27r) [TB]. We therefore obtain the 
spectral determinant for 7 = 0, 

5Dir(0) = det(-A) = exp ( - Cf)i,(0)) = 2L , (7) 

which agrees with the formulation of the spectral determinant of the Laplace operator 
in [m [201 [32]. 



3. 2. Schrddinger operator on a finite interval 

To compute the spectral determinant of the operator —-^ + V{x) for x G [0, L] with 
Dirichlet boundary conditions let us introduce two useful basis for the solutions of the 
differential equation 

(^-^ + V{x)^^{x) = k''ij{x) . (8) 

The first basis shown to be useful in [HI [9] is: /(a;; —k"^) is a solution of (jHl) such that 
/(O; — /c^) = 1 and f{L\ —k"^) = 0, with the second linearly independent solution denoted 
f{x; —k"^) where x = L — x and / satisfies /(O; —k^) = 1 and f{L; —k"^) = 0. For 
a symmetric potential, V{x) = V{L — x) and we have f{x; —k"^) = f{x; —k^). An 
alternative basis was employed in [M]: Uk{x] —k"^) and u{x; —k"^) are solutions satisfying 
u(0; —k"^) = and u'{0] —k'^) = 1, where prime denotes derivation with respect to the 

f A vector potential Ah{xi,) that depends on the position on the bond can be made constant on each 
bond by a gauge transformation. 
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spatial coordinate x, i.e. first argument of the functions. The relation between the two 
basis is given by 

A secular equation whose solutions k are square roots of the eigenvalues of the 
Schrodinger operator is then 

n(L;-P) = 0; (10) 

note that the secular equation may be written in different forms, thanks to the relations 
u(L;-P) = u{L;-e) = -l/f'{L;-e) = -1/ f'{L; -k^), that hold for an arbitrary 
potential. Note, however, that u{x] 7) and u{x] 7) do not coincide in general (apart for 
a symmetric potential V{x) = V{L — x)), nevertheless u{L;'y) = u^L-jj) always holds. 
This may be shown by computing the Wronskian of the two solutions at the ends of the 
interval. 

Example 1: For a free particle where V{x) = 0, 

f{x; —k^) = f{x; —k"^) = — — a^]^^^ ^(^3;- — fc^'j _ ^(•j,. — /^^^ _ gjj^ /j^, _ 

sin kL k 

The Dirichlet spectrum is then given by solutions of the secular equation 

F{k) =u{L- -e) = J]^^^ =lsmkL = 0. (12) 



Example 2: In the case of a linear potential V{x) = co x we find 

TT 

a;~/3 



-Bi{-uj~^/^e)Ai{uj^/^{x-e/uj))] , (13) 



u(x; -e) = ^ \BiU^%L - e/tu)) AiU^^x - kVuj)) 

-M{uj^/'\L - e/u)) Bi(a;^/3(x - z^Vw))] , (14) 

with the Airy functions Ai(z) and Bi{z) |T]. The Dirichlet spectrum is given by solving 
the secular equation 

F{k) = Ai(-w-2/3A;2)Bi(a;i/3(^ _ ^2/^^^ _ Bi{-iu-^/^e)Ai{uj^/^L - k^u)) = . (15) 

Example 3: Another interesting example is the case of a potential of the form V{x) = 
(j){x)^ + <f)'{x). The Hamiltonian 

^5 = -^ + 0(x)' + 0'(x) (16) 



describes so-called supersymmetric quantum mechanics [M]. In this case, exploiting the 
factorization Hs = Q^Q with Q = —-^ + 4>{^)^ it is possible to construct explicitly two 
linearly independent solutions of the differential equation Hsip = 0. We may choose 
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ipoix) = exp (j){y) dy and ^/'i(x) = ^po{x) v^c^fp- obtain the two useful solutions 
for k = 0: 

u{x;0) = ijjo{0)ijjo{x) and u{x; 0) = ipoiL)ipoix) (17) 

We now explain the general method followed in [18] and in the present article. 
A secular equation F{k) = can be related to the zeta-function using the argument 
principle, 

CdUs, 7) = ^ ly + 7)-^^ log (Fiz)) dz , (18) 

where the contour wraps around the positive real axis enclosing the roots of F{z), the 
solutions of the secular equation, and avoiding any poles. Figure [U^a). 




b) 



^ o o-^< o- 



o zero of F( z) 
X pole of F(z) 



-0 ^ — (x) — > 



o zero of F(z) 
X pole of F(z) 



Figure 1. The contours used to evaluate the graph zeta function, (a) before, and (b) 
after, the contour deformation. The two branch cuts are represented by hashed lines. 



In the present case the function F{z) = u{L; —z"^) does not present poles. However, 
as there is some freedom in the choice of F{z), in the following section it is convenient 
to choose a function F for the general graph that will possess poles. Transforming the 
contour C to the imaginary axis z = it, Figure [D^b), we obtain 

Cd.(3,7) = ^ r {t'-^r^^ log {u{L-t')) dt, (19) 
TT at 

where, assuming -F(O) 7^ 0, the segment between 1^7 and — 1^/7 does not contribute to 
the integral as F{z) = F{—z) (which follows from the symmetry of the secular equation); 
comments about the case F{0) = follow in Section [6l Analysis of the integral shows 
that this representation of the zeta-function is valid in the strip 1/2 < JJs < 1. The 
restriction to JJs > 1/2 comes from the asymptotic behavior of u{x] t^) for t — )■ 00. The 
asymptotics of the function u{x; t^) for t — )■ 00 were evaluated in [2], 

\ogu{L-t^) ~ tL-log2t + 0(t°) . (20) 

An explicit example of such an expansion is provided in the appendix, see Eq. (lA.9p . 
To obtain a representation valid for s — )■ we subtract and add the first two terms in 
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the asymptotic expansion. Then upon integration we find, 

Cb-As, 7) = ^ / it' - 7)-' -Ti [log HL- t')) -tL + log(2t)] dt 

TT J at 



vahd in the strip —1/2 < 5Rs < 1. Differentiating 

CDir(0, 7) = - log {u{L; 7)) - log(2v/^) + 1 log 7 = - log (2m(L; 7)) . (22) 

Consequently the zeta-regularized spectral determinant of a wire with Dirichlet 
boundary conditions is 



5Dir(7) = exp ( - CL{0, 7)) = 2n(L; 7) = . (23) 

Example 1: In the absence of a potential we obtain 

, ^ 2 sinh ^fyL , ^ 

^Dir 7 = 7^ . 24 



As 7 — )■ this expression agrees with ((Tj). 

Example 2: Using ( IT3|) . the spectral determinant of H = —-^ + ojx therefore is 
5Dir(7) = [Ai(a;-2/^) B\{uj^/\L + ^/uj)) - Ji\{uj-^'^^) M{u"\L + 7/u;))] • (25) 

Example 3: Using the expression f|T71) we can obtain an explicit formula for the spectral 
determinant of the supersymmetric Hamiltonian for vanishing spectral parameter: 

Myy 



I.e. 

det (^-^ + 02 + 0'^ =2e^o1>i-)'^- j\-^Io^iy)<iydx. (27) 

Note that this explicit formula furnishes another way to obtain the determinant 
det (7 — ^ + V) in terms of one solution among a family of solutions of a Riccati 
equation 4>'{x) = 7 — 4>{xY + V{x). For example, in the free case, V{x) = 0, the Riccati 
problem is solved by 0(x) = ^^7 tanh y/y{x + Xq): we can check that ( 127|) leads to (124|) . 
Formula ( l27l) may also be illustrated in the case (f){x) = ux (harmonic oscillator on a 
finite interval), 

det + w'^' + 00^ = y^e^"^' erf(v^L) , (28) 

where erf(x) is the error function jTB]. The limit L — 00 of fl28p is singular and does 
not approach either a/tt/cj or y^7rj2uj, the determinant of the operator defined on M or 
]R+ respectively. This is not unexpected as eigenf unctions on the interval [0, L] will not 
be square integrable in the limit L — )■ 00. Eigenfunctions of high enough energy always 
feel the presence of the boundary conditions at both ends of the interval and thus differ 
significantly from the infinite space situation. 
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3.3. Dirichlet determinant of a graph 

We are now also in a position to consider the case of a graph with Dirichlet boundary 
conditions at every vertex. The Dirichlet conditions describe a graph where all the 
bonds are decoupled and consequently the spectrum is simply the union of the Dirichlet 
spectra of the individual intervals. Since the zeta-f unctions of each wire are additive 
and the determinants are multiplicative in this case we recover a result of [32] 

To compare with the free case, we use (fTT]) to get S'Dir(7) = 2-^7"^/^ 11^=1 ^inh ^7-^6. 
Taking the limit 7 — we see ( l29i) reduces to the explicit formula S'Dir(O) = 2^ Y[b=i 
which agrees with ([7]). 



4. Secular equation of a Schrodinger operator with general matching 
conditions 

For a general graph let fi, be the solution of the boundary value problem on the interval 

+ h{x,- -P) = eMxb, -k^) , (30) 

such that /fe(0; — A;^) = 1 and fb{Lb', —k^) = 0. fiixj;] —k"^) will denote the complimentary 
solution on the same bond where the function vanishes at the initial vertex of h and is 
unity at the terminal vertex, i.e. /^^(O; —k"^) = 1 and fj;{Lb; —k"^) = o||| fb and are a 
pair of linearly independent functions on the interval satisfying the eigenvalue problem. 
Hence the component of the wavefunction of the Schrodinger operator (jlj) with energy 
k"^ on bond b can then be written 

M^b, -k') = Cb fbixb, -e) e'""^^^ + c-J-^ix-,; -k^) e^^^^^ . (31) 

Matching conditions at the graph vertices are specified by a pair of 2B x 2B matrices 
A and B via 

AV' + BV' = , (32) 

where 

V' = (^i(O), . . . , ^b{0).MLi). . . . , MLB)f , (33) 
= (Di^i(O), . . . , DbMO), DriPiiL,), . . . , DT^B{LB)f . (34) 
Db := — lAb is the covariant derivative and in ( 133|) we consider ^b a function of Xb so 
Xfe = at o{b) and Xb = Lb at t{b). ij) is therefore the vector of inward pointing covariant 
derivatives at the ends of the intervals. 

The following theorem of Kostrykin and Schrader which classifies all matching 
conditions of self-adjoint realizations of the Laplace operator [23] also applies to the 
Schrodinger operator. 

X Notice that it is convenient to consider /j- as a function of xj- = — x^. 
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Theorem 2. The Laplace operator with matching conditions specified by A and B is 
self-adjoint if and only if (A, B) has maximal rank and AM^ = MA^ . 

See [28] for an alternative unique classification scheme for general vertex matching 
conditions. Note that the general matching conditions we employ can be approximated 
by ornamenting the vertices of the graph with subgraphs whose vertices have 6-tjpe 
matching conditions [6]. 

We consider only local matching conditions where the matrices A and B relate 
values of functions and their derivatives on the intervals where they meet at a vertex 
and where the matrices are independent of the metric structure of the graph, namely 
the bond lengths and the spectral parameter k. The restriction on the form of matching 
conditions is necessary in order to employ the argument principle later. 

Representing the components of an energy eigenfunction on G following (13T]) we see 
that = (ci, . . . , cs, cy, . . . , c^)"^ and V' = M{-k'^)il) where M^-k"^) is the 2B x 2B 
matrix 

Ma, = Sa,t /^(O; -e) - S,-, fL{U- -fc2)e'^^^^ . (35) 

The vector ?/? of coefficients of the two linearly independent solutions on each bond 
characterizes the wavefunction. Substituting in (132!) we see that for xj) to define a 
wavefunction of energy fc^ satisfying the matching conditions we require 

det (A + BM(-P)) = , (36) 

which defines a secular equation for the graph Schrodinger operator. 

Example 1: For a free particle Vb{xb) = and we can write fb explicitly, 

,2N smk{Lb-Xb) 2n -k cos k{Lb - Xb) , . 

fbixb] -k = — , so f^ixb, -k = — . 37 

sm kLb sm kLb 

Consequently the matrix M takes the form 

where cot kL = diagjcot kLi, . . . , cot kLs} and esc kL is defined similarly. 

5. Zeta function of a Schrodinger operator with general matching 
conditions 

Let us define the function 

F{z) := det (A + MM{-z'^)) (39) 

so the secular equation ( 136|1 reads F{k) = 0. As previously the zeta function is written 
as a contour integral 

C(^,7) = ;^ [ {z' + ^r^^ log F (z) dz , (40) 
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where the contour C encloses the zeros of F and avoids its poles, Figure [D^a). Making a 
contour transformation to the imaginary axis, Figure [U^b), the integral breaks into two 
parts 

C(s,7) = Cim(s,7) + Cp(5,7) , (41) 
where we denote by the integral along the imaginary axis and (p is the series that 
results from subtracting the residue {zq + •y)'^ at each pole zq of F. Writing z = it the 
imaginary axis integral becomes 

Cim(3,7) = ^^ r{t'-i)-'^\ogF{it)dt, (42) 

TT dt 

which converges in the strip 1/2 < < 1. 

As poles of F lie at eigenvalues of the graph Schrodinger operator with Dirichlet 
boundary conditions we see that Cp(s,7) = CDir(s,7), which we evaluated previously. It 
remains to formulate the imaginary axis integral so that the zeta function given in (l42l) 
is well defined to the left of the line = 1/2. 

5.1. Large t asymptotics of F{it) 

To extend the zeta function to the left of the line 3f?s = 1/2 we require the large t 
asymptotics of F{it) = det(A + BM(t^)). The t — )■ oo asymptotics of fl are analysed in 
the appendix where we find 

oo 

/^(0;t2)^_i + ^3^_^.(0)r^-^ (43) 

i=i 

and fl{L}y]t^) vanishes exponentially. The functions S6j(xb) are determined iteratively 
from the potential Vb{xh) via a recursion relation (1A.4I) . 

Using we can write the expansion of M{t'^) as t — )■ oo up to exponentially 
suppressed terms. 

oo 



Hence 



where 



F(it) ~ det (^A + B(-tI + D(t))) (45) 



D{t) = 5^t"'diag {si,,(0), . . . , SB,,(0), . . . , s^/Q)} . (46) 

i=i 

The expansion of F{it) fits with the free particle case investigated in [18] where, in the 
absence of a potential, D{t) = and F{it) ~ det(A — tM). 

Using fH5|) and the recursion relation flA.4p the expansion of F{it) can be carried 
out to arbitrary order using a computer algebra package. In general the expansion will 
take the form 

oo 

F{it) ~ det (a + B(-tI + D{t))^ = J2 Cj , (47) 

i=o 



Spectral determinants on graphs 



11 



where Cq = det B. We denote by cn the first nonzero coefficient in the expansion of F{it), 
so = if det B 7^ 0. To evaluate the zeta-regularized spectral determinant we need a 
representation of the zeta function valid about s = 0. To find this we only require the 
leading order in the expansion. Subtracting and adding the leading order behavior we 
obtain an expression for the imaginary axis integral valid in the strip — l/2<3fJs<l. 

Cim(^,7) = / (t^-7)-^;TTlog(F(it)t^-2>^)dt + -(2i?-iV)7"^ (48) 

Correspondingly 

Ci'm(0,7) = log - log (F(iy^)7^) - ^^y^ log 7 

= logC7v-logdet(A + BM(7)) . (49) 



6. Zeta-regularized spectral determinant 

Collecting the results so far we can now formulate the spectral determinant, 5'(7) = 
exp ( — CDir(05 7) ~ Cim(05 7))' of ^ graph with general local vertex matching conditions, 
we obtain, 

S(^) = det (A + BM(7)) where S^Ul) = TT 777T^ • (^0) 

This is the statement of Theorem [TJ The 7-dependent part of this expression coincides 
with Desbois' result [9] obtained by a different approach. If we choose Dirichlet matching 
conditions at the graph vertices then A = I2B and B is a matrix of zeros, see ( 132|) . 
Hence 023 = 1, which is the only non-zero term in the large t expansion of F{it), and 
det (A + BM(7)) = 1, consequently f l50|l reduces to fl29|) as required. 

F{z) vanishing at zero. If -F(O) = the contour integral will pick up contributions 
associated with the behavior at z = 0. While -F(O) does not vanish generically, when 
it does vanish it is straightforward to compensate for this by adjusting the function 
F. Assume F[z) ~ z"^^ as 2; — )■ for P G We may define a new function 

F[z) = F{z)/z'^^. Then F{z) = is a secular equation but F does not vanish at zero. 
Evaluating the zeta function using F the procedure is identical and we find, 

5(7) = det (A + BM(7)) . (51) 

Cn 7 

The same technique was used, for example, to evaluate the zeta function of the Laplacian 
on a star graph with Neumann like matching conditions in [18]. To determine P in a 
particular case one requires the asymptotic expansion of /^(O, t^) and fl,{Lb, t^) as t — > 0. 
We now discuss the spectral determinant in two particular cases. 



f The even power is due to the symmetry of F. 
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Functions continuous at the vertices. We consider local vertex matching conditions 
for which the wavefunction is continuous at the vertices, namely the 5-type conditions. 
At a vertex v the local 5-type matching conditions relate the functions and covariant 
derivatives at the ends of the intervals meeting at v. One explicit choice of rUy x ruy 
matrices to encode such matching conditions is 



(52) 



where A„ is the strength of the coupling at v. We note that det(At,— tB„) = —{Xy + niyt). 
By relabeling the 2B ends of the intervals we may write the matrices A and B describing 
matching conditions on the whole graph so the matrices Ay and B^, for v = 1, . . . ,V 
appear as blocks on the diagonal. We recall that the factor cn appearing in the 
spectral determinant is the coefficient of the highest power of t in the expansion of 
det(A + B(— tl + -D(t)), where D{t) is a diagonal matrix where each non-zero element is 
a series in inverse powers of t. Hence we see that for a graph with 6-tjpe conditions at 
all the vertices and any potentials on the bonds 





/ 


— A„ 





■ 


■ M 




/ 1 


1 ■ ■ 


■ ^ M 






-1 


1 


■ 


■ 







■■ 


■ 


V 







-1 


1 ■ 


• 


and B,„ = 





■■ 


• •■■ 




v 








■ 


■ l) 




V 


■■ 


■ ■■■ 0/ 



det(A + BM(7)) ~ det(A 



l{-tI + D{t)) ~ 

t—>oo 



n 



-rriy 



t^ + 0{t 



Thus we find = ^ and cy = Htjl^"^'')- Consequently 

^ ^ det(A + BM(7)) 



n 



-rur. 



(53) 



(54) 



reproducing the result in [52] . Note that the choice of a wavefunction continuous at the 
vertices allows one to introduce vertex variables and simplify the spectral determinant by 
replacing det(A + BM(7)) by a matrix coupling vertices of smaller dimension [21 [HI I 



Functions whose derivative is continuous at the vertices. Another interesting case of 
matching conditions, studied in [3T| 132] . are the so called 5^-type matching conditions 
corresponding to a wavefunction whose derivative is continuous at the vertices [12]. A 
possible choice of matching matrices at a vertex v is 

/ 1 1 ■ ■ ■ 1 1 \ / -fly ■■■ 0\ 






V 












and 



-1 





So det(- 



0/ \ 

-{fly + jniy), therefore 



1 
-1 1 








1 / 



det(A + BM(7)) t^B-^JJl 



-rrir. 



t—^oo 



for /it, = V f . 



(55) 



(56) 
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The exponent = 2B — V and the coefficient C2b-v = Y[v(~^v)- again obtain 



Comparison with the conjecture of W3^ . The zeta-regularized spectral determinant was 
obtained in |32] in the two particular cases considered above with a different method. 
Based on the observation that both sets of matching conditions led to Eq. (154|) it was 
conjectured that this result is of greater generality (Eq. (61) of this reference). In light 
of the present work we now inject two observations. 

• The matrices A and B describing boundary conditions are not unique. Any 
transformation (A, B) — )■ ([/A, UM) where f/ is a matrix such that detU ^ 0, 
leaves the nature of the matching conditions invariant. Therefore the general 
expression of the zeta-regularized spectral determinant should be invariant under 
such a transformation. It is clear that this is not the case for the form (!54ll 
conjectured in [32] whereas our main result Theorem [1] is invariant. 

• The conjecture was based on a continuity argument with respect to variations of 
A and B. However, in general, cat may depend on the particular choice of the 
matching conditions rather than simply the topology of the graph. In particular, it 
is not possible to rule out, a priori, that for a given choice of matching conditions 
and potentials a term in the series expansion (H71) vanishes. We have seen in the 
previous pair of examples that the exponent N does indeed vary with the nature 
of the matching conditions. 
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Appendix. Large t asymptotics of /'(a;; t^) 
f{x;t^) is the solution of 



on the interval [0,L] with f{0]t'^) = 1 and f{L\t^) = 0. To find the asymptotics 
of /'(x;t^) in the limit t — oo we let 0t(x) be a general solution of (lA.ip and call 
St{x) = -plog0t(x). St satisfies the differential equation 




(A.l) 



S[{x)=t' + V{x)-S^{x). 
Consequently, in the limit t ^ oo, St has an expansion of the form 



(A.2) 



oo 




(A.3) 
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The functions Sj{x) are determined by the recurrence relation 

Sj+i{x) = I s'j{x) + Sk{x)sj-k{x) I (A.4) 

\ fc=0 



where S-i{x) = ±1, so{x) = and si{x) = ±^V{x). The two alternative signs in the 
expansion produce exponentially growing and decaying solutions of flA.ll) . If we denote 
with Sf{x) the functions with s_i(a;) = ±1 respectively then we may write the large t 
behavior of f{x\t^) as a linear combination of the exponentially growing and decaying 
parts, 

f{x] t^) = exp / du + A~ exp / (u) du . (A.5) 

Jo Jo 

Imposing the boundary conditions we require A'^ = 1 — A~^ and as it turns out 
, =F exp Sr (u) du 

A = L L • (A-6) 

exp S^^ (u) du — exp S^. (u) du 
Consequently differentiating ( ]A.5|) we find, 

nx;t')= (A.7) 

(x) exp {J^ S;- (u) du + (u) du) - jx) exp{J^ ju) du + (u) du) 

exp <S^{u) du — exp S^^ (u) du 

We are interested in the t — t- oo asymptotics of both /'(0;t^) and From 
(1A.7P we see that at x = L, /'(x; t^) vanishes exponentially. At x = 0, up to 
exponentially suppressed terms, we have the asymptotic relation 

oo 

f'{0-t')^Sr{0) = -t + J2s,{0)t-^ . (A.8) 

Example: In the case of the linear potential V{x) = ux considered earlier the first 
eight functions Sj in the expansion of constructed iteratively from the recurrence 
relation are shown bellow. 

s_i(x) = -1 So{x) = Si{x) = ~x 

S2{x) = -— Ssix) = —X S^{X) = —X 



Consequently 



CO 
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